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ABSTRACT 

We present an analysis of the secular variability of the longitudinal magnetic field in 
the roAp star 7 Equ (HD 201601). Measurements of the stellar magnetic field Be were 
mostly compiled from the literature, and append also our 33 new B^. measurements 
which were obtained with the 1-m optical telescope of Special Astrophysical Obser- 
vatory (Russia). All the available data cover the time period of 58 years, and include 
both phases of the maximum and minimum B^. We determined that the period of the 
long-term magnetic Be variations equals 91.1 ± 3.6 years, with i?e(max) = -|-577± 31 
G and Be(min) = -1101 ± 31 G. 
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1 INTRODUCTION 

The Ap star 7 Equ (HD 201601, BS 8097) is one of the 
brightest objects of this class, with the apparent luminosity 
V = 4.66 mag. The exact spectral type of this object is A9p 
(SrCrEu subclass). The magnetic field of 7 Equ has been 
studied for more than 50 years, starting from October 1946 
(see Babcock 1958). The longitudinal magnetic field Be of 
this star does not exhibit periodic variations in time scales 
typical of stellar rotation, 0.5 — 30 days. Such a variability 
of the Be field was observed in most Ap stars. The above ef- 
fect is commonly interpreted as the result of stellar rotation 
(oblique dipole model). 

The first measurements by Babcock (1958) showed that 
the value of the longitudinal magnetic field Be of 7 Equ 
was positive in 1946-52, and approached nine hundred G. 
From that time on the value of Be slowly decreased and even 
changed sign in 1970/71. One could interpret the magnetic 
behavior of 7 Equ either as secular variations, or variations 
caused by extremely slow rotation. If the latter picture is 
correct, then the corresponding magnetic and rotational pe- 
riods are in the range from 72 to 110 years (Bonsack & 
Pilachowski 1974; Leroy et al. 1994; Bychkov & Shtol' 1997; 
Scholz et al. 1997). 

The behavior of the Be field in 7 Equ was investigated 
by many authors in the second half of the twenty th century. 
For this research we compiled Be observations published by 
Bonsack & Pilachowski (1974), Scholz (1975; 1979), Borra 
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& Landstreet (1980), Zverko et al. (1989), Mathys (1991), 
Bychkov et al. (1991), Bychkov & Shtol' (1997), Scholz et al. 
(1997), Mathys & Hubrig (1997), Hildebrandt et al. (2000), 
Leone & Kurtz (2003) and Hubrig et al. (2004). 

We included in this paper our unpublished magnetic Be 
measurements which were obtained during the past seven 
years. All the new magnetic observations showed, that the 
slow decrease of the Be field in 7 Equ apparently reached the 
minimum in 1996-2002 and has actually started to increase. 

In this paper we determined the accurate parameters 
of secular variability of 7 Equ: the period Pmag, the am- 
plitude and the time of zero phase for Be variations, which 
were approximated by a sine wave. We support the hypoth- 
esis that the long-term Be variation in 7 Equ is a periodic 
feature. Possible origin of this variation cannot be uniquely 
determined, see discussion in Section [Kjof this paper. 



2 OBSERVATIONS AND DATA PROCESSING 

We have performed spectropolarimetric observations of Zee- 
man line splitting for 7 Equ at the Coude focus of the 
1-m optical telescope (Special Astrophysical Observatory, 
Russian Academy of Sciences). Zeeman spectra were ob- 
tained with the echelle spectrograph GECS (Musaev 1996). 
We have put the achromatic analyser of circularly polarised 
light in front of the spectrometer slit. Images of the Zeeman 
echelle spectra were recorded from CCD detectors in stan- 
dard FITS format. Final reduction of the archived spectra 
was performed with the standard MIDAS software (Monin 
1999). 
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Effects of instrumental polarisation on Be measure- 
ments obtained with this instrument were investigated by 
Bychkov et al. (1998, 2000). 

Table Q presents the full set of our Be measurements 
of 7 Equ (total 33 Be points). The meaning of the first 3 
columns is obvious. The fourth column gives the number 
of spectral lines which were used for the measurement of Be 
for a given exposure. Time length At of the exposure (in 
min) is given in the last column of Table Q 

On average, the value of a single Be number listed in 
Table was obtained after averaging of Be measurements 
obtained in 500-1300 spectral lines. Standard deviation as^ 
for the resulting value of Be was computed in the standard 
manner as the error of an arithmetic mean value. 

Errors ct_b^ determined in the above way reached rather 
low values in several observations listed in Table Q In 
2005/2006 we plan to verify the reality of such cr_B^ by a spe- 
cial program of Be observations. Actually we accept these 
errors bona Sde and note the following properties of our Be 
measurements. 

The referee pointed out that a few pairs of Be measure- 
ments of one night in Table ^difi'er by only a few G, which 
is substantially less than the corresponding standard devia- 
tion (Ts^ . We can explain this only as a purely random effect, 
and do not see any reason for it either in the acquisition of 
observational data or their reduction. 

Secondly, series of measurements taken within a few 
nights generally show a scatter of the order of 100 G, which 
is much higher than the standard errors cr_B^ in Table0 The 
latter are of the order of 20 — 30 G, and such a discrepancy 
suggests that our standard deviations are systematically un- 
derestimated, and are in fact of the order of 100 G. On the 
other hand, such a scatter of ~ 100 G is not inconsistent 
with the short-term variability of light and the longitudinal 
magnetic field Be in 7 Equ in time scales of minutes or above 
it. 

Leone & Kurtz (2003) recently discovered periodic vari- 
ations of the longitudinal magnetic field Be in 7 Equ over 
the pulsation period of this star, Ppuis = 12.1 min. The es- 
timated amplitude ABe = 240 G for this period, therefore, 
these variations at least can contribute to the scatter of our 
Be points collected in Table 

Study of the rapid periodic Be variations on a time 
scale of minutes was also presented in Bychkov et al. (2005b) 
for 7 Equ. They did not found conclusive evidence of such 
variations above the noise level at ~ 240 G. 

We also performed spectral analysis of the full set of 
298 Be time series from years 1946-2004. We concluded that 
there are no short-period field variations with periods above 
ca. 1 day, but were not able to extend our analysis for shorter 
periods, see Section 4 of this paper. 



3 MAGNETIC PERIOD OF 7 EQU 

Magnetic observations presented in Table represent com- 
pletely new data. They cover time span of ca. 7 years and 
include the phase when the effective magnetic field Be in 7 
Equ apparently reached its minimum value, and then the 
slow decrease of Be observed in the recent « 50 years has 
been reversed. This fact is of extraordinary importance, be- 
cause it allows one for a fairly accurate determination of the 



Table 1. Measurements of Be in 7 Equ (HD 201601). 



JD 2400000. -f 


Be (G) 


o"s„ (G) 


N 


At (min) 


49648.323 


-1045 


21 


706 


30 


49648.345 


-1315 


26 


755 


30 


49649.229 


-1463 


37 


576 


30 


49649.257 


-1159 


31 


656 


30 


49932.424 


-1317 


26 


691 


60 


49932.469 


-1317 


26 


675 


60 


49933.460 


-1316 


26 


700 


60 


49933.507 


-1317 


29 


704 


60 


50023.158 


-1291 


22 


501 


40 


50023.189 


-1380 


23 


650 


40 


50066.128 


-1539 


26 


718 


40 


50066.157 


-1611 


62 


532 


40 


51533.1229 


-1014 


16 


966 


30 


51533.1451 


-1011 


14 


701 


30 


51535.1847 


- 902 


16 


955 


40 


51535.2153 


- 901 


19 


855 


40 


51536.1069 


- 670 


18 


821 


30 


51536.1285 


- 642 


24 


508 


30 


51888.166 


-1069 


18 


847 


30 


51888.190 


-1092 


20 


1353 


30 


51889.103 


- 890 


20 


847 


30 


51889.126 


- 865 


20 


817 


30 


51890.142 


- 742 


21 


770 


30 


52163.3000 


- 845 


19 


833 


30 


52163.3201 


- 855 


19 


732 


30 


52164.2861 


- 956 


16 


947 


30 


52164.3076 


- 967 


16 


914 


30 


52165.2812 


-1061 


17 


835 


40 


52165.3111 


-1029 


16 


991 


40 


52186.2229 


- 922 


17 


1085 


30 


52186.2451 


- 942 


17 


1055 


30 


52187.2673 


- 882 


16 


1072 


30 


52188.2395 


- 908 


18 


838 


30 



magnetic period and the amplitude of Be variations in 7 
Equ. 

We have compiled the set of 298 observations of the Be 
field in 7 Equ, scattered in the literature, and appended our 
measurements. These data cover the time period 1946-2004 
(58 years). They are displayed in Fig. Note, that the Be 
measurements obtained by Babcock (1958) apparently cover 
the phase of the maximum longitudinal magnetic field in 7 
Equ. 

The set of Be measurements analysed in this paper is 
rather heterogeneous. The data have been obtained by sev- 
eral different observers over a long time period using various 
instruments and techniques, and it is impossible to estimate 
or test credibly their systematic and random errors, partic- 
ularly for the earliest observations of the longitudinal mag- 
netic field in 7 Equ. 

Therefore, we arbitrarily assumed that systematic er- 
rors of the Be observations are equal to zero. In other words, 
all the Be points for 7 which were found in the literature are 
fully compatible. 

Random errors of individual Be points frequently were 
given in the source papers, and are denoted by vertical bars 
in Fig.Q These errors were not directly available for the ear- 
liest photographic measurements by H.W. Babcock (1958) 
and Bonsack & Pilachowski (1974). We adopted here the 
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estimated error for Babcock's data equal 238 G, and 151 G 
for Bonsack & Pilachowski. These numbers were obtained in 
our thorough reanalysis of the earliest papers dealing with 
measurements of stellar magnetic fields, of. Section 3.1 in 
Bychkov et al. (2003). 

Determination of the period and other parameters of 
the apparent magnetic variability for 7 Equ was performed 
in the following manner. Assuming that the run of the ob- 
served longitudinal field Be with time T can be approxi- 
mated by a sine wave 



Be(T) ^Bo + Bi sin 



27r(r-To) _ TT 
P 2 



(1) 



we determined all four parameters: the period P, the average 
field Bo, the amplitude Bi and the time of zero phase To 
using the iterative technique of nonlinear fitting. 

Starting values of P, Bq, Bi, Tq and their standard 
deviations were found by our computer code for the non- 
linear least squares method (Bychkov et al. 2003). The fi- 
nal values and their errors were then computed with the 
public domain code "nlfit.f", which is designed for curve 
and surface fitting with the Levenberg-Marquardt procedure 
(ODRPACK V. 2.01 subroutines). The code is available at 
the site www.netlib.org. 

Fitting of a sine wave to all the 298 Be points with er- 
rors as in Fig. gave very poor results with the for ^ 
single degree of freedom /'^ = 18.0420. Such fits are un- 
acceptable, and in case of 7 Equ the poor fit is the result 
of underestimated errors of many Be points. Many Be ob- 
servations presented in Fig. have very low errors, which 
sometimes are less than 20 G. Our new Be points, which are 
collected in Tabled also are of such a high formal accuracy. 

We cannot judge, whether an apparent scatter of Be 
points in Fig. is due to unrealistic error estimates or the 
intrinsic short-term variability of the longitudinal magnetic 
field in 7 Equ. The estimated random error of Be points 
about the starting sine wave equals to 213 G. For the final 
fitting of a sine we assumed that all the 298 points have 
identical errors of 213 G. 

Final values of the fitted parameters and their standard 
deviations a for the sine phase curve are given below. 

Pmag = 33278 ± 1327 days = 91.1 ± 3.6 years 

To = JD 2417795.0 ± 1057. 

Bo = - 262 ± 22.4 G 

Bi = -f 839 ±22.1 G 

r = -0.524 ±0.043 

In other words, a parameter range from —a to +a is just 
the true 68% confidence interval for this parameter. 

The above fit of a sine wave with uniform errors of 213 
G is very good, with x^/^^ = 1.0134. The effect of inhomo- 
geneity in the Be time series plus the possible existence of 
rapid magnetic variability in 7 Equ were compensated by 
the increase of the random error, and neither should infiu- 
ence the above parameters of secular magnetic variability in 
7 Equ. 

The standard parameter r was defined for the oblique 
rotator model of an Ap star. It is related to the angle (3 
between the magnetic dipole axis and the rotational axis, 
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Figure 1. The longitudinal magnetic field Be for 7 Equ in years 
1946-2004. 



and the angle i between the rotational axis and the line of 
sight (Preston 1967): 



cos /3 cos i — sin /3 sin j _Be(min) 
cos f3 cos i + sin /3 sin i Be (max) 



(2) 



Parameters Be (min) and Be (max) of the Be sine wave for 7 
Equ are given by 



Be (max) =Bo±Bi 
Be (min) = Bo - Bi 



± 577 ±31.4 G 
-1101 ±31.4 G 



Note, that the meaning of Be (max) and Be (min) for 
use in Eq.|21is different: Be (max) denotes there the value of 
magnetic intensity which has the higher absolute value, and 
Be (min) has the lower absolute value. In this way we ob- 
tained the value of r for 7 Equ equal to r = 577/ ( — 1101) = 
-0.524. 

Bychkov et al. (2005a) presented an extensive catalog of 
the magnetic phase curves and their parameters for 136 stars 
on the main sequence and above it. We quoted there the pre- 
viously estimated period for 7 Equ, Pmag = 27027"^, which 
was obtained on the basis of a shorter series of Be data. 
This paper and the new, more accurate Pmag = 33278^* rep- 
resents a major revision of the previously known magnetic 
period of 7 Equ. 



4 SEARCH FOR ADDITIONAL MAGNETIC 
PERIODS IN 7 EQU 

Significant scatter of the observed points in the long-term 
run of Be(T) in Fig. suggests the search for short-term 
periodicities. We applied the strategy of prewhitening to the 
set of available Be measurements, and removed the principal 
sine-wave variations from the data. Prewhitened data were 
then analysed with the method developed by Kurtz (1985), 
and with his Fortran code (Kurtz 2004). 

Such a search for peaks in the Be amplitude spectrum 
of 7 Equ in this paper was restricted to trial periods higher 
than 1 day. This is because many of the earlier magnetic 
observations for this star either have poorly determined the 
time of measurement, or have long times of exposure (see 
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Figure 2. Amplitude spectrum of tlie _Be time series for 7 Equ, 
years 1946-2004. 



e.g. Babcock 1958). The star 7 Equ exhibits rapid nonra- 
dial pulsations and the corresponding with the period 
Pmag ~ 12.1 miu (Leone & Kurtz 2003) and, possibly, with 
simultaneous shorter periods (Bychkov et al. 2005b). None 
of them were analysed in this paper. 

We have identified two additional periods of statistically 
low significance in the range Pmag > I'', see Fig.|5| 

Pi = 348.07 days, amplitude = 122 G 
P2 = 23.44 days, amplitude = 110 G 

Both peaks in the amplitude spectrum in Fig. |5| exhibit low 
signal to noise ratio, with noise level at ca. 80 G. The period 
Pi is close to 1 year. Since most of the existing Pg observa- 
tions for 7 Equ were performed in months July-November, 
then the peak Pi in the amplitude spectrum represents a 
false period which most likely reflects the average 1-year 
repetition time in the acquisition of the existing magnetic 
measurements. 

We believe that the peak P2 in the amplitude spectrum 
of the Pe field of 7 Equ is the random effect of a pure noise. 
The peak is very narrow, in fact, it only appears in a single 
bin of a very dense discrete frequency mesh. 

Kurtz (1983) discussed the possible existence of the pe- 
riod of « 38 days in his photometric observations of 7 Equ 
in 1981. That period was of low probability, but possibly 
could be identified with the real rotational period in this 
star. We do not confirm the existence of the 38 day period 
in long-term Pe observations of 7 Equ, see Fig. |5| 



5 DISCUSSION 

There exist three possible explanations for the observed 
long-term behavior of the longitudinal magnetic field in 7 
Equ: 

1. Precession of the rotational axis (Lehmann 1987). 

2. Solar- like magnetic cycle (Krause & Scholz 1981), 

3. Rotation with the period of 91.2 years. 



The Ap star 7 — HD 201601 in fact is a binary sys- 
tem. One can assume, that the gravitational force from the 
secondary companion can cause precession of the Ap star. 
As the result, the angle between the rotational axis and the 
direction towards the Earth varies periodically. Therefore, 
changes of the aspect can in principle cause apparent varia- 
tions of the longitudinal magnetic field Pe or the amplitude 
of its variations. 

Effects of precession in long-period Ap stars were stud- 
ied by Lehmann (1987), who showed that the oblateness of 
stars caused by the rotational or magnetic flattening is not 
adequate to produce observable precession effects. The only 
exception was 52 Her, where the observed behavior of the 
star could be interpreted as a precessional motion. 

The above considerations indicate that the precession 
theory does not convincingly explain Pe variations in this 
star. 

The idea by Krause & Scholz (1981) that we actually 
observe the solar-like magnetic cycle in 7 Equ in which the 
global magnetic field reverses its polarity, cannot be easily 
verified by the existing observations of the global longitudi- 
nal magnetic field Pe. Moreover, one can note that such an 
idea requires the existence of a mechanism in the interior of 
7 Equ which ensures the transfer of huge magnetic energy 
into electric currents and vice versa. Note that the required 
efficiency of such a mechanism and the amplitude of mag- 
netic field variations in 7 Equ is ca. four orders of magnitude 
larger than that in the Sun in a similar timescale. 

Following the widely accepted picture of an Ap star, we 
believe that the magnetic field of 7 Equ can be approximated 
by a dipole located in the center of the star. The dipole is 
inclined to the rotational axis of 7 Equ. We assume that the 
magnetic field is stable and remains frozen in the interior 
of a rotating star at the time of observations, i.e. during at 
least of 58 years. Therefore, slow variations of the Pe field 
in 7 Equ are caused by an extremely slow rotation, in which 
case our Pmag = Prot = 33278''. Such an explanation is 
supported to some extent by polarimetric measurements by 
Leroy et al. (1994). 

We plan to perform high accuracy polarimetric mea- 
surements of 7 Equ with the new version of MINIPOL. The 
device was constructed to measure the angles and the degree 
of linear polarisation of stellar radiation, and will be oper- 
ational at the Special Astrophysical Observatory in 2006. 
We also expect that shall be able to verify the extremely 
slow rotation of 7 Equ measuring the rate of change for the 
polarisation angle of stellar radiation. 



6 SUMMARY 

The Ap star 7 Equ (HD 201601) exhibited slow and system- 
atic decrease of the longitudinal magnetic field Pe starting 
from 1946, when the global magnetic field of this star was 
discovered (Babcock 1958). We have compiled the full set 
of 298 existing Pe measurements, which consists of the Pe 
data published in the literature and our observations ob- 
tained during recent 7 years. The latter magnetic data (33 
Pe points) were measured with the echelle spectrograph in 
the Coude focus of the 1-m telescope at the Special Astro- 
physical Observatory. Our newest observations showed that 
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the longitudinal magnetic field of 7 Equ reached its local 
minimum and started to rise in 1998-2004. 

All the available data cover the time period of 58 years 
(1946-2004) and include both phases of the maximum and 
minimum Be. Assuming that the secular variability of the 
Be field is a periodic feature, we determined parameters of 
the magnetic field curve in 7 Equ and give the value of its 
period, P — 91.1±3.6 years, with the zero phase (maximum 
of Be) at To = JD 2417795.0 ± 1057. Sine-wave fit to the Be 
phase curve yields Be (max) = +577 ± 31 G and Be(min) = 
-1101 ±31 G. 

Spectral analysis of the 58-year long Be time series es- 
sentially do not show the existence of shorter periods, down 
to trial periods of ~ 1 day. More specifically, there arc no 
real shorter periods in the run of the longitudinal magnetic 
field Be with amplitudes exceeding the noise level of 80 G. 
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